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■ Abstract 

(N . 

. We calculate the cross section for the production of two jets with large transverse momenta 

^^O ' in DIS diffractive dissociation for both transverse and longitudinally polarized photons. The 

■ scale which defines the hardness of the Pomcron is found to be . We present analytic 
i-^^l expressions and discuss numerical results relevant for the diffractive dissociation at HERA. 

iu: 
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•1—1 , 

^ ' 1.) The diffractive dissociation of the photon in deep inelastic ep-scattering at HERA has recently 

' attracted much interest. As for the inclusive diffractive cross section attention has been given 

to the question whether the Pomeron which describes the energy dependence of this process is 
closer to the soft Pomeron which describes the high energy behavior of hadron-hadron scattering 
or whether it resembles the hard Pomeron observed in the proton structure function at small x. 
Whereas earlier analysis [^, § of HERA data seem to favor the first alternative, a recent analysis 
IQ] presents evidence for a strong admixture of the hard Pomeron. A natural way to gain further 
clarification is a more detailed analysis of the diffractive final state. Electroproduction of vector 
mesons and photoproduction of J/^' have been analysed Q and compared Q with both the soft 
and the hard Pomeron, and they seem to favor a harder Pomeron. A new class of final states which 
allow to enhance the hard Pomeron component are jets with large transverse momenta ^, The 
simplest configuration are two-jet events with the jets originating from qq pairs; for larger invariant 
masses M one expects events with additional production of gluon jets to become dominant. 
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In this paper we consider the simpler case qq jet production, restricting ourselves to the kine- 
matic region where the hard Pomeron can be expected to dominate (fig. 1). Starting from the DIS 
diffractive cross section formula 
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we calculate the differential cross section for the subprocess ^* + p 
consider both the transverse and longitudinal cross section. 



-|- p for very small t and 



2.) We assume the total energy s = (p + q)'^ to be much larger than the photon virtuality = —q^ 
as well as the invariant mass of the quark anti-quark pair = {q + xjpp)"^, i.e. the Pomeron mo- 
mentum is much smaller than the proton momentum (xjp <C 1). Instead of one can also 
introduce the variable /3 which is defined as /3 = Q"^ /{M"^ + Q^)- For our analytic calculation we 
set the momentum transfer t to zero, since the cross section strongly peaks at t = 0. In order to 
compare with data we will add the t-dependence later on by hand. 

The new element in our approach is the fixed transverse momentum k of the (anti-)quark which 
enters as a second hard scale in our calculation and allows to treat the Pomeron perturbatively 
(hard Pomeron). For large s (small xjp) the amplitude of the process is dominated by two-gluon 
exchange, and following the k-factorization theorem |^] we express the amplitude with leading- 
log(l/xjp) accuracy through the unintegrated gluon distribution of the proton. The two gluons 
have to couple to the outgoing quarks in all possible combinations in order to preserve gauge 
invar iance. 

The quark phase space is usually parameterized in terms of the transverse momentum k and the 
longitudinal momentum fraction a according to the Sudakov decomposition: k + xjpp = aq' + 
k^/ {as)p + k (the anti-quark momentum consequently has the parametrization: q — k= (1 — a)q' + 
k^/[(l — a)s]p — k) with q' = q + xp. The invariant mass can now easily be expressed in terms of 
a and k^: = 2{q — k) • {k + xpp) = k^/[a(l — a)]. In the quark anti-quark CMS one finds 
for a = (1 — cos(0))/2 where 9 denotes the angle between the quark and the proton. From the 
expression k^ = a{l — a)M'^ it becomes clear that keeping the mass M fixed (of the order of Q) 
and taking k^ to be small also means that a and 6 become small. As will be seen below the cross 
section is dominated by small k^, i.e. the jets appear most frequently at small angles or phrased in 



different words the jets are 'aligned' (see Aligned Jet Model ||1(_ 

The cross section for the qq production has the following form (for details see refs. [0, 
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where represents the unintegrated gluon distribution of the proton 

J dl''TGixF,f) = xipG{x]p,Q^) . (6) 

This identification of the Pomeron with the gluon structure function, however, has to be taken 
with care. Strictly speaking the function which appears in eq. (^) is not exactly the same gluon 
structure function as in DIS. Namely the longitudinal component of the two gluon lines in fig. 1, 
even at t = 0, have a typical value of the order of xjp and their difference is exactly xjp. In 
the leading-log(l/a;jp) approximation where we cannot distinguish log(l/xjp) and log(2/xjp) the 
function G from (^) has the same xjp-dependence as the gluon structure function. But in the 
HERA-regime we cannot exclude a difference in the absolute normalization. In the following we 
shall identify our function G with the gluon structure function, but in a more ambitious analysis 
one could think of calculating the nonforward gluon structure function using the generalized gluon 
splitting function of p6[| . 

Note the extra in the denominator of the longitudinal cross section in (3) which reflects 
the fact that this is a higher twist contribution. The term [...J/l^ in the integrand of eq. @ is 
roughly constant for 1^ smaller than k'^{M'^ + Q'^)/M'^ and falls off hke 1/1^ for larger 1^. The 
unintegrated gluon structure function J^g also falls like 1/1^ (modified by logarithms) so that 
the dominant contribution in the integration comes from the region 1^ < ]<.'^{M'^ + Q'^)/M'^. We 
approximate the integrand here by taking the limit 1^ = for and integrate over J^g with 

k2(M2 + Q'^)/M'^ as the upper limit. The latter results in X]pG{x]p,k^{M'^ + Q^)/M^), i.e. we 
have extracted the leading log{li^{M'^ + Q'^)/M'^) contribution. The caveat, however, is that this 
approximation underestimates the region of large M^. This becomes clear by first taking the limit 
— > oo and then performing the integration over 1^. In this limit [...] equals 2^(1^ ~ k^), i.e. 
there is no contribution at low 1^. The integrand at large P is roughly falling like 1/1^ and the 
integral is determined by its lower limit k'^ which approximately leads to J-G{x]p,k.'^)/k-^ , i.e. the 
final result also contains the unintegrated structure function. In the next section we present the 
correction for any Af- value which unfortunately cannot be derived from simple arguments as given 
here. Physically the scale h'^{M'^ + Q'^)/M'^ corresponds to the virtuality of the softer of the two 
quarks, into which the photon dissociates. With the approximation described above one obtains: 
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where bt and 6/ are functions of ]VP and and will be given in (11) and (12). Let us discuss 
the main properties of this result. It is seen that the transverse momentum of the outgoing quark 
pair sets the scale of the gluon structure function. This means that the effective slope of the hard 
Pomeron, which determines the rise of the cross section at small xjp should increase with increas- 
ing transverse momentum. The effect is even strengthened by the additional factor {Q'^ + M'^)/M'^ 
which enters the scale and is of the order of two for the Q^,M2-range considered here. From the 
known behavior of the gluon structure function we can furthermore conclude that the power rise 
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of the cross section as a function of is damped for small, but still in a range where the 
perturbative approach is legitimate. Note that (7) and (8) cannot be extrapolated down to k^ = 0. 
The fall-off with increasing is identical for both polarizations. A remarkable property of the 
first term of the longitudinal cross section in (8) is its zero for = M^. A numerical evaluation 
of the above formulas will be performed in part 4.). 



3.) Thus far we have discussed the cross section in the double-leading-log approximation. One can 
in principle do better because the k- factorization formulas eq. (|5|) are more general and allow 
for an evaluation without the leading-log(Q2)-approximation. This is equivalent to using the BFKL 
Pomeron |17] for the unintegrated gluon structure function. We will show that the characteristic 
momentum scale which was identified in the previous section emerges in this model without any 
approximation and that we essentially reproduce the formulas given there, provided this scale is 
large enough. 

The use of the BFKL Pomeron requires a nonperturbative input distribution which cannot 
be determined theoretically. Using only the most general parameterization we find for the BFKL 
Pomeron the following expression 
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with Ag being a nonperturbative scale (Ag — IGeV^), (j)^!^) a slowly varying function of v and 
■y^^u) = ^^^^[2ip{l) - 'ip{l/2 + iv) - ■0(1/2 - iv)] the characteristic eigenfunction of the BFKL 
pomeron. Now we insert the above expression in our basic equations (|5|) and perform the 
^-integration exactly. This leads to 
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where F is the hypergeometric function. The remaining z^-integral can be safely evaluated nu- 
merically but here we restrict ourselves to the discussion of some limiting cases. One observes 
immediately that the scale k2(Q2 _|_ m'^)/M'^ is the relevant parameter, which determines the 
location of the saddle point of the integral. As long as NcOs/tt • logl/xjp is much larger than 
A = log[k2(Q2 _^ m2)/(m2Ao)] we are in the BFKL-hmit and the saddle point lies at 1/ ~ 0. In 
this case one finds the usual BFKL results with a steep rise of the cross section for small xp. 

If, on the other hand, A is much larger than Inl/xjp, we are in the double logarithmic limit 
and the i^-contour has to be shifted close to z/ = —i/2. The saddle point is located at ivg = 
1/2 — y^Ncas/^TrA) ln(l/xip). Evaluating the z/-integral we obtain 
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If we compare these results with those obtained in the previous section, we see that an improvement 
has been achieved in the prefactor of the unintegrated structure function. It is more accurate in 
the transverse case and we have found a correction for the longitudinal part, too, where the zero 
for Q2 = jvf2 disappears. The leading term, however, which contains the structure function itself 
is the same as before. We should stress here again that the scale of the gluon structure function 
k2(Q2 _j_ m2)/m2 already appeared in the single logarithmic approximation and turned out to be 
the most important parameter in this process. 

Finally, let us comment on the region of small k2. The expressions given in (10) and (11) can 
be continued down to k2 = and the integral gives a finite answer. For the transverse cross section 
a detailed saddle point analysis shows that the main contribution comes from the small k2-region 
and thus agrees with the Aligned Jet Model. For the longitudinal cross section on the other hand, 
the additional factor k2 makes the small k2-region less dominant. Hence, for the k2-integrated 
transverse cross section the BFKL model gets a too big nonperturbative contribution and thus is 
not reliable. The longitudinal cross section has less nonperturbative contributions. 

The situation becomes better, if k2 is constrained to be large. This was shown for diffractive 
jet production in a numerical study contained in ||l8[| . From this we draw the conclusion that our 
above analysis, for both cases, is on a solid theoretical basis if only A, which is the logarithm of 
the relevant scale here, is large. Hence we regard the eqs. (|l2|), ([T3| ) as a firm prediction, up to a 
normalization uncertainty, in a range, where k2 is at least of the order of IGeV'^ and where the rise 
of the gluon density as a function of x is visible, the latter condition ensuring that we are in the 
small-x limit. 



4.) Based upon the formulae (|2|), (^), (|l^ ) and (^) we have performed a few numerical cal- 
culations which may apply to HERA data. Throughout our calculations we have integrated over t 
by taking our cross section expression evaluated at t = and multiplying with the formfactor which 
was given by Donnachie and Landshoff |1£]. Furthermore we used the running strong coupling con- 
stant with the scale being given by the scale of the gluon structure function. First we have looked 
into the most striking feature of the cross section formula, the rise of the gluon structure function 
at small xjp. Fig. 2a and 2b show the xjp-dependence of the transverse and the longitudinal cross 
sections, keeping f3 and fixed at 2/3 and 50GeV^, respectively, and integrating over k2 > 2GeV^. 
We have considered different parameterizations for the gluon structure function: the GRV leading 
order and next-to-leading-order parameterizations [^, and the hybrid model of one of us (M.W.) 
|14|. For comparison we also show a prediction of the soft Pomeron M. As expected, the hard 
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Pomeron rises much stronger at small xjp than the soft Pomeron. The upper three curves show 
the overall uncertainty in both xjp-shape and overall normalization of the perturbative approach. 
In order to find the "optimal" version we have calculated the inclusive structure function F2, by 
coupling the gluon to the quark loop and doing the same type of approximations as for the qq 
cross section described in the previous section (here we have assumed that in the kinematic region 
which we are considering the gluon structure function dominates). Figs. 3a-d show the results of 
this computation in comparison with data |21], and we feel that the next-to-leading order GRV 
parameterization gives the best description. After having in this way "selected" which parameteri- 
zation for the gluon structure function is the best one for our purposes, we return to eqs. (12) and 
( [l3| ) and calculate the xjp-dependence, varying now the momentum scale li^{Q^ + M'^)/M'^ which 
determines the "hardness' of the Pomeron structure function. As it can be seen from fig. 2c, there 
is a clear variation in the slope (numerical values for the slopes are given in the figure), but in the 
HERA range this variation is not very strong. So the main experimental signal, for the beginning, 
may be the observation of a rise in eqs. dH) and (H) with a power of 1/xp around 0.75. 

Next we have used our formulae for estimating the ep-integrated cross section for the production 
of two jets with large transverse momenta. We define two cuts; (a) k^ > 2GeV^ and (b) k^ > 5GeV'^, 
and we then integrate over all other variables in the region 

Q2 > lOGeV^, xip< 10-2 and 50GeV < 
W < 220GeV {W'^ = {p + q)'^)- The results for the tr ansverse and the longitudinal cross sections 
are given in Table 1. Typically the longitudinal cross section is of the order 10% of the transverse 
cross section. The other conclusion that may be drawn from this table concerns the dependence 
upon the k^ cut: moving from the cut (a) to cut (b), one loses about 75% of the cross section. 
In fig. 4 we show the dependence on k^, keeping j3 and fixed. On the whole one notices that 





kg = 2GeV^ 




kg = 5GeV^ 






GRV(LO) 


GRV(NLO) 


GRV(LO) 
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31 
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Table 1: Results for total ep-cross sections (in pbarn) of diffractive dijet production for two different 
parameterizations of the gluon density and two different cuts on the transverse momentum of the 
jets. 

the longitudinal cross section (figs. 4 b,d) decreases less rapidly than the transverse one (figs. 4 
a,c). According to our formulae, the transverse cross section falls approximately as 1/k^, whereas 
the longitudinal one decreases more like l/k'^. This simple power behavior is, however, slightly 
modified by the k^-dependence of the gluon structure function (numerical values for the slopes 
are given in the figures). In the region of smaller k^ the cross sections start to become somewhat 
flatter. This is in agreement with what we have already discussed after eq. (P). The most striking 
variation is seen in fig. 4b: in the double leading log approximation the longitudinal cross section 
has a zero at = M^: due to the correction term in ([l3|), this zero is somewhat shifted towards 
smaller (3. What we see in fig. 4b is the remainder of this zero. 

Another quantity of interest is the /3-spectrum for jets with k^ larger than 2Ge\f^ (fig. 5). As 
to the transverse cross section, in comparison with the inclusive diffractive cross for qq production 
which shows the (1 — /3)-dependence we now see a maximum which at not too large lies below 
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/3 = 0.5. In the longitudinal cross section which again is smaller than the transverse cross section 
by a factor of 10 we see the dip which has moved from the value (3 = 1/2 (uncorrected cross section) 
to /? = 0.4 or less. 

5.) In this paper we have collected evidence that the observation of two jets with high trans- 
verse momentum may provide a valuable way of analysing, in the DIS Diffractive Dissociation, 
the hard Pomeron. After having established that the scale for the hardness of the Pomeron is not 
simply but the combination k^(Q^ + M^)/M^ we have illustrated several signals of this hard 
Pomeron. In our numerical estimates we have found that it is not enough to consider that part of 
the phase space where the transverse momenta are strongly ordered. The main signal of the hard 
Pomeron is, most likely, the stronger rise in 1/xjp. But there are also other, more subtle details 
which are typical for the hard Pomeron. A more careful comparison with the predictions of the 
soft Pomeron seems very important. 

Production of two jets clearly represents the simplest case of jet production in DIS Diffractive 
Scattering. In the next step one has to compute higher order corrections to two jet production 
and to generalize to the production of of extra gluon jets which are expected to become important 
especially in the large M-region. Another point of interest is the t-dependence. Whereas the soft 
Pomeron predicts the typical shrinkage for the t-slope, the hard Pomeron is expected, at least for t 
of the order of IGeV^, to have a much weaker energy dependence of the t-slope. Both these issues 
have to be analysed in more detail in order to reach a more complete understanding of the hard 
Pomeron in Diffractive DIS. 
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figs. 2 a,b. We gratefully acknowledge the help of C.Ewerz in comparing our F2 results with the 
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Figure Captions 



Fig. H : A contribution to the amplitude for the process 7* + p — > qq + p. The shaded blob 
represents the unintegrated gluon structure function. For the full amplitude the gluons have 
to be coupled to the quarks in all possible ways. 

Fig. 2a : The Xip-dependence of the cross section da'^ ^ /d(3 for transverse photons. The curves 
are for = 50GeV^, (5 = 2/3 and integrated from 2GeV^ up to the kinematical limit. 
Results are shown for the soft Pomeron exchange model of 0, the perturbative hybrid model 



of 1 14 1 and the model presented in this paper with the GRV leading order and next-to-leading 



order gluon distribution 
Fig. ^ : The same as fig. 2a but for longitudinal photons. 



Fig. 2c : The j;jp-dependence of the cross section • da"^ ^/dp for transverse photons us- 
ing eq. (p!^) and the GRV(NLO)-parameterization for three different kinematical situations: 
Q2 = 80GeV2, p = 2/3, integrated from 2GeV2 to 4GeV2 (solid hue), = SOGeV^, 
/? = 2/3, k2 integrated from 4GeV2 to SGeV^ (dotted line) and = sOGeV^, /3 = 5/6, k^ 
integrated from 2GeV^ to 4GeV^ (dashed line). The numbers in the figure give the slope 
of each curve. The normalization is the integral of the cross section over the xjp-range 
displayed. 

Figs. 3a-3^ : Comparison of our results for the inclusive structure function F2, calculated with the 
same type of approximation as described in section 3.) and two different parameterizations 



of the gluon density [gO[ with data from HI , ZEUS and E665 [21] for different values of Q 



Fig. 4a : The k^-dependence of the differential 7* p-cross section for transverse photons according 
to eq. (0). Values for the kinematical parameters are xjp = 5 ■ 10~^, (3 = 2/3, and is 
varied between 15GeV^ and 45GeV^. The quantity S gives the effective slope of the curves, 
obtained from a numerical fit to a power behaviour ~ (k^)"''. 

Fig. ^ : The same as in fig. 4a but for longitudinal photons (eq. (p!^)). 

Fig. ^ : The same as in fig. 4a but for fixed = lOGeV^ and different values of /3. 



Fig. ^ : The same as in fig. ^ but for longitudinal photons. 



Fig. 5a : The /3-dependence of the cross section da^ ^ /dfi for transverse photons. The curves are 
for xjp = 5 • 10"'^ and three different values of Q^. 



Fig. pq : The same as in fig. 5a but for longitudinal photons 
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